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We present results from studies of high temperature QCD with two flavors of Kogut-Susskind quarks on 16^ x 8 
lattices at a quark mass of amq = 0.00625 and on 24^ x 12 lattices at quark masses airiq = 0.008 and 0.016. 
The value of the croRsovor temperature is consistent with that obtained on coarser lattices and/or at larger quark 
masses. Results are presented for the chiral order parameter and for the baryon number susceptibility. 



In this note we report on two studies of high 
temperature QCD with Kogut-Susskind quarks. 
The HTMCGC collaboration has recently com- 
pleted a study on 16^ x 8 lattices at a quark mass 
of aniq = 0.00625 where a is the lattice spac- 
ing[l]. This project is a follow-on to earlier work 
on the same sized a lattices at am, = 0.0125[2]. 
The MILC Collaboration is in the process of car- 
rying out calculations on 24^ x 12 lattices with 
quark masses aniq = 0.008 and 0.016. The ob- 
ject of these projects is to push the study of high 
temperature QCD towards the continuum limit 
and physical quark masses in order to determine 
the nature of the transition between ordinary 
hadronic matter and the quark-gluon plasma, 
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and to make more reliable calculations of the 
transition temperature and other thermodynamic 
quantities. 

In the Nt = 8 project we have performed sim- 
ulations for gauge couphngs in the range 5.45 < 
6/(7^ < 6.00. A total of 1,000 trajectories (molec- 
ular dynamics time units) were generated for each 
value of the coupling studied except = 5.475 
for which 1,450 were obtained. In Figure 1 we 
show the average value of the chiral order param- 
eter, (V'V'), and the baryon number susceptibil- 
ity, X, obtained from these runs as a function of 
As expected, x is small in the low temper- 
ature regime since a baryon must be created or 
destroyed to change the baryon number, while it 
is large in the high temperature regime since there 
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Figure 1: a) The chiral order parameter and b) 
the baryon number susceptibihty as a function of 
6/ff2 on 16^ X 8 lattices at am, = 0.00625. 



it is only necessary to create or destroy quarks to 
change the baryon number. Just as on coarser 
lattices [3], the baryon number susceptibility pro- 
vides an excellent signal for the crossover. We 
estimate the crossover value of the coupling at 
this quark mass to be 6/<jf^ = 5.487(13) 

For Nt = 12 we are studying gauge couplings 
in the range 5.65 < 6/5^ < 5.85 for quark masses 
airig = 0.008 and 0.016. The autocorrelation 
times are large, and we plan to obtain approx- 
imately 2,000 trajectories at each coupling by 
the time this project is completed. At present 
we have over 1,000 trajectories at all couplings 
at which we are running except = 5.85, 

arrig = 0.008. However, the results presented be- 
low must be treated as preliminary. 

In Figure 2 we plot the chiral order param- 
eter, (V'V')i a-s a function of for both of 
the masses which we are studying, arug = 0.008 
and 0.016. ^AV' also show linear extrapolations to 
zero quark mass obtained from pairs of points at 
arrig = 0.008 and 0.016. For some couplings ex- 
trapolations are shown from more than one pair 
of points. The bursts are extrapolations from 
two points with the same coupling, but different 
masses, while the fancy crosses are extrapolations 
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Figure 2: The chiral order parameter, (V'V') a-s a 
function of 6/g^ on 24^ x 12 lattices. 



from points with different couphngs. Using re- 
sults for the chiral order parameter and for the 
Polyakov loops, we estimate that for arrig = 0.008 
the crossover occurs at = 5.71(2). It is pre- 
mature to estimate the position of the crossover 
for arUg = 0.016. 

In Figure 3 we plot Tc/rrip as a function of 
{rriT^Imp)'^. Here, is the temperature at the 
pseudocritical point determined in simulations 
run at finite quark mass: 



Tjmp = l/(Ntamp) 



(1) 



In order to produce this figure we have made fits 
of m„ and nip as a function of gauge coupling and 
quark mass, from spectrum calculations with two 
flavors of Kogut-Susskind quarks. These flts al- 
low one to extrapolate am„ and airip to the quark 
mass and pseudocritical value of the gauge cou- 
pling of the thermodynamics simulations. The 
non-Goldstone pion is used in this plot to point 
out the large flavor symmetry breaking. The 
line on the left hand side of the plot marks the 
physical value of (mjr/mp)^. Results with Wil- 
son quarks are included for comparison. The 
curved line running through the Nt = & Wflson 
point shows the path it might fall on while stay- 
ing within error bars. The Nt = 12 point is for 
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arriq = 0.008. It is quite preliminary, and it would 
be premature to attribute any significance to the 
fact that it is slightly larger than values of Tc/nip 
at larger lattice spacings. Rather, we consider it 
quite striking that there is such good agreement 
between the Kogut-Susskind results form Nt = 4 
to Nt = 12, and between them and the Wilson 
results for Nt = 6 and 8. 
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Figure 3: Tc/rrip as a function of (mjr/m)p)^- The 
non-Goldstone pion is used for Kogut-Susskind 
quarks. 

The mass spectrum fits can also be used to cal- 
culate the baryon number susceptibility in phys- 
ical units. This quantity enters into the QCD 
equations of state. It also plays an important 
role in models of nucleosynthesis, hadronization 
and heavy ion collisions. In Figure 4 we show 
the present world's data. The insensitivity of the 
results to lattice spacing is encouraging. 

The scaling properties of thermodynamic quan- 
tities are important for determining the nature of 
the transition to the quark-gluon plasma. The 
results presented here arc consistent with the ex- 
pectation that for two fiavors of quarks, there is a 
second order phase transition for zero quark mass, 
and a rapid crossover for finite quark masses [4]. 
This topic is discussed in some detail in DeTar's 
review talk[5]. 

The HTMCGC Cohaboration's calculations at 
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Figure 4: The baryon number susceptibility in 
physical units as a function of temperature. 



Nt = 8 were carried out on the Connection 
Machine CM2 at the Pittsburgh Supercomputer 
Center. The MILC Cohaboration's Nt = 12 cal- 
culations are being performed on the Connec- 
tion Machine CM5 at the National Center for 
Supercomputer Apphcations, Intel Paragons at 
the San Diego Supercomputer Center and San- 
dia National Laboratory, and Cray T3D's at the 
Pittsburgh Supercomputer Center and the Na- 
tional Energy Research Supercomputer Center. 
The work was supported by the NSF and DOE. 
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